Role of the first coordination shell in determining the equilibrium structure and 

dynamics of simple liquids 
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The traditional view that the physical properties ol a simple liquid are determined primarily by 
its repulsive forces was recently challenged by Berthier and Tarjus, who showed that in some cases 
ignoring the attractions leads to large errors in the dynamics [L. Berthier and G. Tarjus, Phys. 
Rev. Lett. 103, 170601 (2009); J. Chem. Phys. 134 214503 (2011)]. We present simulations of 
the standard Lennard- Jones liquid at several condensed-fluid state points, including a fairly low 
density state and a very high density state, as well as simulations of the Kob-Andersen binary 
Lennard- Jones mixture at several temperatures. By varying the range of the forces, results for 
the thermodynamics, dynamics, and structure show that the determining factor for getting the 
correct statics and dynamics is not whether or not the attractive forces per se are included in the 
simulations. What matters is whether or not interactions are included from all particles within the 
first coordination shell (FCS) - the attractive forces can thus be ignored, but only at extremely 
high densities. The recognition of the importance of a local shell in condensed fluids goes back to 
van der Waals; our results confirm this idea and thereby the basic picture of the old hole- and cell 
theories for simple condensed fluids. 



I. INTRODUCTION 

Ever since van der Waals [l| in 1873 formulated his 
theory "the continuity of the gaseous and liquid states" , 
the general understanding of simple fluids has been that 
repulsive and attractive forces give rise to separate mod- 
ifications of the ideal-gas equation of state and to sepa- 
rate contributions to the free energy. The harsh, repul- 
sive forces reduce the available volume and determine the 
structure of the fluid and thus its configurational entropy; 
the weaker, longer-range attractive forces give rise to an 
energy effect that reduces pressure and energy compared 
to those of an ideal gas at the same temperature and 
density. With respect to the effect of the longer-range 
attractive forces van der Waals imagined that he could 
(quoting from Ref. 0) "... define an clement of volume 
in a liquid which is small compared ..." with "... the 
range of the intermolecular force, but large enough for 
it to contain sufficient molecules for us to assume that 
there is within it a uniform distribution of molecules of 
number density p" . He assumed that the main effect of 
the attractive forces originates from molecules within this 
local sphere and estimated this sphere to have a range of 
3-6 A|2]. Below we present evidence that this local vol- 
ume, although it does not have a uniform distribution of 
molecules, may be identified with the first coordination 
shell (FCS). 

It is well known that the van der Waals equation of 
state reproduces the qualitative behavior of the fluid 
state [3- 5]. Van der Waals' idea about the role of the 
repulsive and attractive forces lies behind Zwanzig's high- 
temperature expansion [6], in which the contribution to 
the free energy from the long-range forces is expressed 
in powers of the inverse temperature, the reference high- 
temperature system being a system with infinitely strong, 
purely repulsive forces. The success of this perturba- 



tion expansion was demonstrated by Longuet-Higgens 
and Widom 0, Barker and Henderson Q, and soon af- 
ter by Weeks, Chandler, and Andersen (WCA) in their 
seminal paper "Role of repulsive forces in determining 
the equilibrium structure of simple liquids" [9(. 

In perturbation theory a fundamental problem is how 
to separate the strong repulsive from the weaker, long- 
range attractions. One possible separation was proposed 
by Barker and Henderson, who marked the separation at 
the distance where the potential is zero. WCA demon- 
strated, however, that by choosing instead the separa- 
tion at the potential minimum one obtains a much bet- 
ter agreement between the particle distributions of sys- 
tem and reference system. Doing so separates the forces 
into purely repulsive and purely attractive forces, and 
one may say that the original idea of van der Waals is 
here captured in its purest form. Numerous refinements 
of perturbation theory have since appeared; of particular 
relevance to the findings reported below are the works of 
Ree et al. from 1985 [l0| and of Hall and Wolynes from 
2008 [ll[ , who studied the effects of a shifted- forces cutoff 
placed at a distance that scales with density in the same 
way as the radius of the FCS. An excellent summary of 
perturbation approaches up until 1976 can be found in 
Barker and Henderson's classic review (l2j . 

The present work is motivated by two recent papers of 
Berthier and Tarjus. In 2009 they showed [l3| that the 
WCA approximation gives much too fast dynamics for 
the Kob-Andersen binary Lennard- Jones (KABLJ) [l4| 
viscous liquid. They concluded that, while the attractive 
forces have little effect on the liquid's structure, these 
forces affect the dynamics in a "highly nontrivial and 
nonperturbative way" . This led Pedersen et al. to in- 
vestigate whether it is possible to reproduce both statics 
and dynamics of the KABLJ liquid by purely repulsive 
potentials different from the WCA potentials [15| ■ This is 
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indeed possible using inverse power-law (IPL) potentials; 
even the thermodynamics is slightly better predicted us- 
ing IPL potentials than using the WCA approximation 
to the true LJ potentials [la ]. Altogether, this indicates 
that it is not the presence of attractive forces per se, 
which is responsible for getting the correct dynamics, 
leaving open the question why the WCA approximation 
fails so dramatically for the KABLJ liquid's viscous dy- 
namics. Very recently the problem was reconsidered by 
Berthier and Tar jus [16j . From a careful numerical inves- 
tigation of the standard LJ and the KABLJ liquids they 
conclude that in the viscous regime the WCA approxima- 
tion gives results that are quantitatively and qualitatively 
different from the true potentials, a finding "which ap- 
pears to contradict the common view that the physics of 
dense liquids is dominated by the steep repulsive forces" . 
Moreover, they find that "a key aspect in explaining the 
differences in the dynamical behavior of the two models 
[WCA vs. true potential] is the truncation of the inter- 
action potential beyond a cutoff" . 

This last point, in conjunction with recent results of 
ours on the shifted- forces cutoff [17j . inspired to the 
present investigation that systematically studies the ef- 
fect of neglecting the forces beyond a shifted- forces cutoff. 
We conclude below that the crucial point is to include all 
interactions from the FCS; the effect of interactions be- 
yond the FCS on the thermodynamics can be taken into 
account to a good approximation by first-order pertur- 
bation theory. 

The radius of the FCS for a fluid of simple spheri- 
cally symmetrical molecules can either be defined as the 
distance where the radial distribution function has its 
first minimum or, following van der Waals 0, as the dis- 
tance from a particle for which the mean density within 
a sphere with this radius equals the overall particle den- 
sity. For the LJ system these two distances are almost 
identical (Fig. [TJ; for the KABLJ system there is a small, 
but insignificant difference. 

A condensed fluid may be thought of as similar to a dis- 
torted crystal with holes. Along this line of thinking, the 
importance of the FCS in condensed fluids was acknowl- 
edged long ago by Frenkel Il8| . Eyring [l9j], and Lennard- 
Jones and Devonshire [2CJ, |2l| . Hole theories and lattice 
theories for liquids are discussed, e.g., in Refs. [H. 
The findings reported below confirm van der Waals' idea 
of the central role of the forces from particles within a 
local shell where the number density equals the overall 
density, as well as the basic ideas behind the hole and 
lattice theories for condensed fluids. 

Section II presents our results for the standard 
Lennard- Jones system, Sec. Ill investigates the highly 
viscous KABLJ mixture. Section IV summarizes briefly 
the consequences of our findings for the general picture 
of simple liquids and for the role of perturbation theory 
for understanding the physics of simple liquids. 
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FIG. 1: Radial distribution and local mean density in a con- 
densed LJ fluid at temperature T — 1.00 and density p = 0.85 
(black straight line). With red <?lj(?")> with green the mean 
density within a sphere with radius r and a particle at its 
center. The minimum of ghj(r) at r — 1.55<r indicated by 
an arrow (corresponding to 5.28A for Ar) coincides with the 
radius for which the local density equals the mean density 
p = 0.85. 



II. SIMULATIONS OF THE LENNARD-JONES 
FLUID WITH VARYING CUTOFFS 

The basic idea behind the standard perturbation ex- 
pansion for fluids is that the particle distribution is de- 
termined primarily by the repulsive part of the potential. 
In this generally accepted picture the effect of the longer- 
range attractive part of the potential can be obtained 
from the first (mean-field) contribution in the expansion. 

For a given state point with density p and temperature 
T the validity of the mean-field approximation may be 
tested by determining structure and dynamics of the fluid 
when the spatial range of the interactions is varied in the 
simulations via the cutoff. We make the model as simple 
as possible. The basic pair potential is the Lennard- Jones 
(LJ) function 

u(r) = M(r/<?r 12 (r/<?)- 6 } ■ (1) 

The range of interactions is limited by introducing a 
shifted-forces (SF) cutoff at a radius denoted by r c 
[l7L l24j . A SF cutoff has the pair force go continuously 
to zero at r c , which is obtained by subtracting a constant 
term as follows 

, ,v jfhj(r) -/Lj(r c ) if r<r c 
JSF{r) = < (2) 
10 it r > r c . 

This corresponds to the following modification of the po- 
tential: Usp(r) = u LJ (r) - (r - r c )u' jrl (rc)- ULj(r c ) for 
r < r c , usf{t) = for r > r c . Reference fL7l showed that 
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a SF cutoff gives more accurate results than the tradi- 
tionally used shifted-potential (SP) cutoff |25| . 

Neglecting the forces beyond r c leads to a difference 
in the thermodynamics, dynamics, and structure of the 
system. After introducing a cutoff the pair distribution 
function g(r) deviates from the "true" ffLj(r). In the 
mean-field approximation, the energy and pressure in- 
clude the effect of the neglected forces as a mean contri- 
bution, calculated from the particle distribution of the 
reference system (which ignores the long-range forces). 
For the pressure, p(p 7 T), the mean-field approximation 
is 1 

2.7T f°° 

p{p,T) = P T + —p 2 g(r)f hJ (r)r 3 dr. (3) 

Any difference between the mean-field approximation of, 
e.g., the pressure and the true value is caused by a differ- 
ent particle distribution originating from the exclusion of 
the forces beyond r c . By varying the cutoff we can thus 
test the validity of the mean-field approximation. The 
validity of the first-order perturbation expansion may be 
investigated by determining the difference between the 
(radial) distributions 

A.g(r) =g(r) -ffu(r). (4) 

The potential energy and the pressure deviate from the 
correct values (denoted by ?7lj and plj) by the terms 

poo 

AU{p,T)/N = 2irp / Ag(r)u hJ (r)r 2 dr, (5) 
Jo 

and 

n p oo 

Ap(p,T) = —p 2 ^ Ag(r)f LJ (r)r 3 dr, (6) 

In molecular dynamics these integrals are obtained di- 
rectly as time and particle averages of the explicit inter- 
action energies, u(rij), and forces, /(r^j), between pairs 
of particles, i and j . In standard perturbation theory to 
first order the Hclmholtz free energy difference between 
the system with a cutoff and the true LJ system is given 
by 

AF(p,T) = AU(p,T). (7) 

The perturbation expansion is asymptotically correct 
at high temperatures for systems with infinitely strong 
purely repulsive forces Q.i 

Five state points of the standard single-component LJ 
condensed fluid were investigated [26]. The five state 
points are: 

1. p = 0.85, T = 1.00 (liquid at medium temperature); 



2. p = 0.85, T — 0.65 (liquid at low temperature and 
almost zero pressure); 

3. p — 2.5, T — 100 (condensed fluid at very high 
density and temperature); 

4. p = 0.60, T = 1.5 (low-density condensed fluid just 
above the critical temperature); 

5. p = 1.05, T = 1.00 (fee solid). 

The fluid systems studied consisted of N = 2000 parti- 
cles, the fee solid system of N = 12 3 = 1728 particles (us- 
ing always periodic boundary conditions). For each value 
of the cutoff the systems were simulated for at least 16 
million time steps (« 160ns in Argon units). Below we 
present data for the fluid state points (l)-(4), for which 
results from simulations using various values of the cutoff 
are compared to results for the true LJ system, which is 
here defined by the cutoff r c = 4.5ct (these data are de- 
noted by LJ). Results for the crystalline state point (5) 
confirm the general physical picture arrived at. 

A. Thermodynamics and dynamics 

Figure Ufa) shows the change in percent of the poten- 
tial energy as a function of the cutoff obtained by Eq. (5) 
at state points (1) (red) and (2) (green), Fig. &b) shows 
the change in pressure using Eq. (6) for the same simu- 
lations. Figure [3Ja) shows the difference in self-diffusion 
constant AD for the LJ fluid at state points (1) (red) and 

(2) (green), Fig. [3Jb) shows the same for state points (1), 

(3) , and (4). The arrows mark the locations of the first 
minimum of the radial distribution function, delimiting 
the FCS. 

For both thermodynamics and dynamics we find the 
following: There is virtually no effect on structure and 
dynamics of ignoring the forces beyond the standard cut- 
off used in most LJ simulations, r c = 2.5a. This is not 
surprising in view of the fact that the interactions are 
very small at such large distances. There is a small effect 
on the dynamics of neglecting the forces from distances 
1.5cr < 2.5cr. Interestingly, the accuracy of the simula- 
tions "recovers" and the simulations generally agree well 
with those of the true LJ system when r c w 1.5cr. This 
is where g^j (r) has its first minimum (marked by arrows 
in Figs. [5] and [3]), delimiting the first coordination shell 
(FCS). As forces are gradually removed for r < 1.5(7, 
however, both thermodynamics and dynamics begin to 
deviate significantly from those of the true LJ system. 
In particular, the self-diffusion constant increases signifi- 
cantly. Note that using the WCA cutoff r c = 2 1 / 6 (at the 
potential minimum), the energy obtained by the mean- 
field approximation deviates only a few percents from the 
correct value. This indicates that even modifications of 
the potential that lead to only quite small changes in free 
energy (Eq. (7)) may significantly affect the dynamics. 

To summarize the findings so far: 1) Simulations of 
the LJ liquid's thermodynamics and dynamics are fairly 
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FIG. 2: Changes of potential energy A U and pressure Ap as 
functions of the cutoff (a shifted-forces cutoff is used through- 
out the paper). In each figure the arrow marks the location 
of the first minimum of g(r), which delimits the first coordi- 
nation shell (FCS) (each figure should have two arrows, but 
they virtually coincide and only one is shown), (a) AU/Ulj 
in percent as a function of the cutoff in the LJ fluid with 
density p — 0.85. The data in red are for T = 1 at which 
U/N — —5.7396 (state point (1)), the data in green are for 
T = 0.65 at which U/N = -6.1242 (state point (2)) (b) 
Ap as a function of the cutoff in the LJ fluid with density 
p — 0.85. The data in red are for T = 1 (state point (1)) 
where plj = 2.0769, the data in green are for T — 0.65 where 
Plj = -0.0912 (state point (2)) g^. 



reliable when the cutoff is placed at <?(r)'s first minimum; 
2) the small energy changes - and thereby small free en- 
ergy changes - induced by removing even a minor part 
of the attractive forces within the FCS lead to a signifi- 
cant increase of the self-diffusion constant. This suggests 
the following. For a condensed, uniform fluid the forces 
on a given particle from particles within its entire FCS 
play crucial roles for the dynamics, whereas the forces 
from particles outside the FCS have little influence and 
their contribution to the thermodynamics may be taken 
into account to a good accuracy via mean-field methods. 
There is a small effect on the thermodynamics and, in 



FIG. 3: Relative changes of the self-diffusion constant AD 
at various cutoffs, (a) AD/ADlj in percent as a function of 
the cutoff in a LJ fluid system at the density p — 0.85. The 
data in red are for T = 1 where Dlj = 0.0497, the data in 
green are for T = 0.65 where Dlj = 0.0258 (the arrow marks 
the location of the first minimum of g(r)). (b) AD/ADlj for 
different densities. The red curve is for state point (1), the 
black curve is for the high-density, high-temperature state 
point (3) (p,T = 2.5, 100). The blue curve is for state point 
(4) (p,T = 0.60, 1.5). The three arrows mark the locations of 
the first minima of g(r) for the respective state points. 



particular, on the self-diffusion constant by also including 
the forces just outside the FCS: by including the attrac- 
tive forces from the interval 1.5a < r < 2. Oct the energy 
decreases by one percent, whereas these attractive forces 
lower the self-diffusion constant by up to four percent for 
state point (2). 

With the van der Waals picture and perturbation the- 
ory in mind the following question now arises: Is it im- 
portant that the forces within the first coordination shell 
includes some attractions? The answer to this question 
is no. Recall that a purely repulsive inverse power-law 
system does exist that has almost the same radial distri- 
bution function as the L J system [28|, [29[ . Since the po- 
tential of the mean force w(r) on a particle is given by the 
radial distribution function ^30], w(r) — — kBTln[g(r)], 
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the mean force on a LJ particle from the particles in its 
FCS does not depend directly on the sign of these forces 
but merely on g(r). 

In order to demonstrate that the dynamics of a LJ 
fluid is obtained correctly even when only purely repul- 
sive parts of the LJ potential lie within the FCS, we sim- 
ulated the LJ system at such a high density that the 
entire FCS is within the "WCA range" delimiting the 
range of the repulsive forces (r < 2 1 / 6 ): p — 2.5, state 
point (2). Even at this state point a cutoff at 3lj(^)'s 
first minimum (r c = 1.09cr) gives good results for self 
diffusion (Fig. [3]Jb)). This shows that inquiring into the 
respective roles of repulsive versus attractive forces is less 
central than focusing on the role of the FCS for simple 
condensed liquids' thermodynamics and dynamics - the 
main conclusion of the present paper. Figure [3Jb) shows 
the relative deviation of the self-diffusion constant D for 
state points (1) and (3), demonstrating that this behavior 
is independent of whether forces are attractive or repul- 
sive. Included in this figure is also the relative change 
in self-diffusion constant for a LJ fluid at state point (4) 
(blue dashes). This low-density state point does not re- 
ally have a well-established coordination shell, and the 
diffusion does not show the same strong dependence of 
the forces within the first shell of nearest neighbour; nev- 
ertheless, this state point shows the same trend as the 
other state points. 



B. Structure 
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Figure Ufa) shows ffLj(r) at the high-density, high- 
temperature state point (3) (black) along with the corre- 
sponding distribution at state point (1) (red). The 3lj(0 
of the high-density state point (3) is almost a simple scal- 
ing of <?lj(?~) of state point (1) to a shorter distance [29| . 
with maximum shifted from r sa 2 1 / 6 er at state point (1) 
to r » 0.75cx at state point (3). As mentioned already, the 
first minimum in <?Lj(f) at state point (1) is shifted from 
r r; 1.5cr to a value at state point (3) of r w 1.09er, i.e., to 
a distance where all the particles within the FCS inter- 
act with the central particles by purely repulsive forces. 
Shown with solid points are values of g(r) for the systems 
where the forces beyond the FCS at state points (1) and 
(3) have been neglected. The two sets of distributions 
(line: ghj(f) and dots: g(r)) are practically identical. 
This demonstrates that forces beyond the FCS play a 
negligible role for the distribution, no matter whether 
these forces are attractive or repulsive. The green dashes 
show the corresponding distributions obtained when only 
the forces for distances shorter than the distance at the 
first maximum of ffu(r) are taken into account in the 
dynamics, i.e., using for state point (1) r c = 2 1 / 6 er and 
for state point (3) r c — 0.75a. 

In perturbation theory the (radial) distribution of par- 
ticles around a target particle is to a good approximation 
given by the distribution obtained from the short-range 
part of the pair interactions 12]. For "WCA" perturba- 



FIG. 4: (a) Radial distribution functions for the two state 
points (1) ((p,T) = (0.85,1.00)) and (3) ((p,T) = (2.5,100)) 
for different cutoffs. The black line gives <7lj(t") at state point 
(3), the red line gives gLj(r) at state point (1). The solid 
points give g(r) at (1) with r c = 1.5<r and at (3) with r c — 
1.09cr, respectively, i.e., for a simulation that includes only 
forces from particles within the FCS. The dashed green curves 
give g{r) for (1) with r c = 2 1//6 <r and (3) with r c = 0.75cr, 
respectively, i.e., keeping only forces up until the maximum 
of g(r). (b) Ag(r) for different values of the cutoff at state 
point (1). With green is shown results for a (WCA) cutoff at 
the minimum of the LJ potential (r c = 2 1 ^ 6 cr). The red line 
is for a cutoff at the first minimum of <;lj = 1.55<r, the blue 
line is for r c = 2.5<r. 



tion theory the unperturbed system neglects all attrac- 
tive interactions. Although the two distributions (green 
dashes and red lines of Fig. 0Ja)) agree reasonably well 
and lead to a fairly good thermodynamic description us- 
ing the first-order perturbation corrections to the ther- 
modynamic functions, there are systematic differences. 
To some extent, by neglecting the attractive forces in the 
dynamics there is a cancellation of errors in the radial dis- 
tribution function in the integrals of Eqs. (H}-©. Figure 
|U[b) highlights Ag(r) at state point (1) for three values of 
r c . With green line is shown the differences for the WCA 
cutoff (purely repulsive forces), the red line gives Ap(r) 
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after including all forces within the FCS (the correspond- 
ing g(r) is shown in the right part of Fig. E^a) with green 
dashes and red point, respectively). For reference, the 
blue line gives Ag(r) for the standard shifted-potential 
cutoff at 2.5a. The figure demonstrates that substantial 
improvement is obtained by including all particles within 
the FCS in the dynamics, although there are still small 
differences. The effect of these small differences cancels, 
however, to a much larger extent than when using the 
purely repulsive WCA reference system, and it leads to 
an excellent overall agreement for both thermodynamics 
and dynamics (Figs. [5] and [3]). 



III. RESULTS FOR THE KOB- ANDERSEN 
BINARY LENNARD-JONES VISCOUS LIQUID 

The previous section showed that what matters for 
simulating with high accuracy the condensed LJ liquid's 
thermodynamics, dynamics, and structure, is to ensure 
that all interactions within the FCS are included in the 
dynamics - independent of whether or not there are any 
attractive forces coming from particles within the FCS. 
We moreover showed that accurate results are obtained 
by placing the cutoff right at <?(r)'s first minimum, which 
provides a convenient and obvious definition of the bor- 
der of the FCS. The 2009 paper by Berthier and Tar- 
jus [l3[, which reopened the discussion of the role of the 
attractive forces in simulations, focused on the highly 
viscous Kob-Andersen binary Lennard- Jones (KABLJ) 
liquid For this system the WCA cutoff leads to a 

dynamics that is up to two orders of magnitude too fast. 
Thus the highly viscous KABLJ system is particularly 
sensitive, and it is therefore important to test whether 
the dominance of the forces coming from particles within 
the FCS identified in the LJ simulations applies also for 
the KABLJ liquid. 

The KABLJ liquid is a mixture of 80% large (A) and 
20% small (B) LJ particles (obb = 0.88ctaa) at the den- 
sity pa AA —l .2 and with a very strong AB attraction [lij ]. 
The A particles dominate the overall dynamics, and the 
(few and small) B particles are to a large extent slaves 
of the structure set by the A particles. For this reason, 
when we refer below to the FCS of the KABLJ system, it 
means always the first coordination shell of the A parti- 
cles delimited by the first minimum of gAA{r). Likewise, 
cutoffs of all interactions are given in units referring to 

(TAA M- 

The mean-square displacement of a particle in a vis- 
cous fluid separates into two regimes: the short-time bal- 
listic regime in which the particle vibrates within its FCS, 
and the long-time diffusive regime that reflects the occa- 
sional escape from the shell. The importance of the range 
of the forces for viscous fluids is investigated below by 
simulating the KABLJ system of N = 1000 particles at 
P a AA = l-^ an d T — 0.45. In order to determine the time 
correlations in this highly viscous fluid state accurately 
the simulations were very long, typically 160 million time 



1000 




time 



FIG. 5: (a) Log-log plot of the mean-square displacement 
(msd) of the A particles in the KABLJ fluid as a function of 
time for different values of r c . With black is shown the "true" 
msd for r c = 4.5<taa; green: r c = 1.12a a ^ [3~H|: light blue: 
r c = 1.2<taa; blue: r c = 1.3<taa; magenta: r c = 1.4<taa; red: 
r c = 1.5<taa; yellow: r c = 2.0<jaa and dark blue: r c = 2.5<jaa 
(there are no visible differences for r c > 1.5<taa). (b) Log- 
log plot of the msd of the B particles in the KABLJ fluid as 
a function of time for different values of r c . With black is 
shown the "true" msd for r c = 4.5<taa; green: r c = 1.12a a ^ 
[3l| : light blue: r c — 1.2<taa; blue: r c = 1.3<taa; magenta: 
r c = 1.4<taa; red: r c = 1.5ctaa; yellow: r c = 2.0<taa an d 
dark blue: r c — 2.5<taa (there are no visible differences for 
r c > 1.5<taa). Dark small dashes: the msd of the A particles 
for r c — 4.5(Jaa, also shown in (a) with black. 



steps. 

The mean-square displacements (msd) for different 
cutoffs are shown as functions of time in Figs. 0(a) and 
(b), which gives msds of the large (A) and the small 
(B) particles, respectively. Both figures show the same 
qualitative behavior. The displacements change drasti- 
cally and systematically by increasing the cutoff from 
r c = \.Ylo a fi to r c = 1.5cr aa, but for r c > 1.5oaa there 
are no visible differences - here both A and B particles 
are trapped for a long time, which is independent of the 
long-range attraction. 
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r(cut) 

FIG. 6: (a) Incoherent intermediate scattering function 
F s {q,t) at the wave vector qoAA = 7.25 for the A particles 
in the KABLJ fluid at T — 0.45 and density pa AA =l.2 as a 
function of time for different values of r c . With black is shown 
F(q,t) for r c — 4.5<taa. Green: r c — 1.12a a ,p [3l]], the line is 
for T = 0.45 and green dashes is for T = 0.3125; light blue: 
r c = 1.2<taa; blue: r c = 1.3<taa; magenta: r c = 1.4a aa; red: 
r c = 1.5ctaa; yellow: r c — 2.0<taa and dark blue: r c = 2.5<taa 
(the differences for r c > 1.5<taa are within the uncertainties), 
(b) The relaxation time r of the KABLJ fluid at T = 0.45. 
With read line and error bars are shown r for the A parti- 
cles obtained from Eq. (8). The self-diffusion constants D^ 1 
(green line) and D^ 1 (blue line) are scaled to agree with r for 
r c = 4.5<taa- 

The A particle self-intermediate scattering function at 
qa aa = 7.25, 



^(?>*) = ^E e *' ( * ,(t) ~ r ' (0)) / > ( 8 ) 

is traditionally used for probing the dynamics of the 
KABLJ viscous liquid. This function is shown in Fig. 
inja) for various cutoffs. The conclusion is the same 
as for the behavior of the mean-square displacements: 
forces from particles beyond r ~ 1.5ctaa play little role 



for F s (7.25,t). The two green curves in Fig. EJa) are 
for a WCA cutoff (r c = 2 1 / 6 ct Qi/3 ) [H. The green line 
is for T = 0.45 and the green dashes (almost invisible) 
is for a lower temperature, T = 0.3125. At this lower 
temperature the self-intermediate scattering function and 
the self-diffusion constant for the KABLJ- WCA system 
agree with those of the KABLJ system at T = 0.45. The 
dynamics of the KABLJ- WCA mixture at T = 0.3125 
is indistinguishable from that of the KABLJ system at 
T = 0.45, which shows that the two systems agree after 
a simple temperature shift. A similar scaling behavior 
has been reported for simple LJ fluids [32j |. (However, it 
should be noted that the KABLJ- WCA system is more 
prone to crystallization because the exothermic binding 
between solute and solvent molecules is not present. It is 
possible to maintain the stability of this viscous mixture 
by keeping the exothermic attraction between solute and 
solvent particles [Hj|.) 

Relaxation times were calculated from F s (7.25,t) via 

/>oc />oo 

t= tF s (7. 25, t)dt/ / F s {7. 25, t)dt. (9) 
Jo Jo 

Results are shown in Fig. [BJb) along with those obtained 
from the self-diffusion via the slopes (6D) of the mean- 
square displacements using t(D) oc a 2 /D, where the re- 
laxation times t(D)a and t(D)b were scaled to agree 
for r c = 4.5(taa- Within the statistical uncertainties the 
three relaxation times agree and show the same depen- 
dence on r c . These findings for the dynamics of the highly 
viscous KABLJ fluid confirm that it is the forces from dis- 
tances below 1.5(7aa that give the correct high viscosity 
and long relaxation time of the viscous state. 

Inspection of the radial distribution functions reveals 
that the FCS, which determines the dynamics, is the one 
established by the A particles. Figure [TJa) shows the 
three radial distribution functions for the KABLJ sys- 
tem. The threshold value r c ks 1.5ctaa for the KABLJ 
fluid is almost the first minimum of the radial distribu- 
tion function <7aa(j) of the A particles, which appears 
at r c = 1. 425cr aa- Note that this is at a distance quite 
different from the first minima of the two other radial dis- 
tribution function; in fact the B particles do not create a 
FCS of their own. 

In a recent paper Rehwald, Rubner, and Heuer demon- 
strated that KABLJ systems of less than one hundred 
particles have the correct (bulk) thermodynamics and 
diffusivity 34]. Tripathy and Schweizer further very re- 
cently showed that in the activated barrier hopping the- 
ory based on naive mode-coupling theory and nonlinear 
Langevin equation, the short-range interactions deter- 
mine the single-particle dynamics and the physical nature 
of the transiently arrested state (in fact, even for non- 
spherical particle fluids) [35[. The present finding that 
it is enough to know the interactions within the FCS fits 
nicely into the results of these two recent works. 

We end this section by briefly digressing to discuss in 
more detail the roles of the A and B particles. That it 
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FIG. 7: (a) Radial distribution functions, g a ,p(r) for the 
KABLJ fluid at T = 0.45. Red line: gAA(r); green: gAs(r) 
and blue: <?bb(V)- 5 ) Radial distribution functions, gBs(r) 
for the KABLJ fluid T = 0.45. Black line: £ BB (r) for 
r c = 4.5<taa; green line: gsB{r) for r c = \.Yla a ^ [jy]]; 
red line: <?bb(?") for r c = 1.5ctaa; magenta line: gBs(r) for 
r c = 2.0<taa (calculated with the cutoff at r c = 4.5ctaa). 



is the FCS around the A particles, which determines the 
properties of the KABLJ fluid, can be deduced from the 
dependence of the radial distribution function of the B 
particles. All three radial distribution functions depend 
on r c , but whereas 3aa( 7 ') and ^abM vary only little 
with r c for r c > 1.5oaa ~ and in a way similar to that 
of the simple LJ-fluid - the distribution of the B parti- 
cles is more sensitive to r c for values larger than 1.5oaa- 
This may seem to contradict the above conclusion, but 
in fact, as we shall now argue, it merely demonstrates 
that the A particle FCS determines the properties of the 
KABLJ fluid. The sensitivity of <?bb(?") to the trunca- 
tion of forces beyond r c = 1.5<jaa is clear from Fig. 
[7J Only when r c > 2<taa is the correct radial distri- 
bution function obtained. Nevertheless, the msd of the 
B particles (Fig. [5](b)) and the corresponding relaxation 
time (blue line in Fig. [Hb)) agree well with the KABLJ 
msd and r already for r c = 1.5<taa- Thus the attrac- 



tive forces between pairs of B particles at distances in 
the interval 1.5<taa < r < 2.0ctaa changes the distri- 
bution <7aa(t)) but have no influence on the dynamics. 
A plausible explanation of this apparent contradiction 
is that a pair of B particles can be bound to the same 
A particle within this particle's FCS and rearrange be- 
tween different local minima within this FCS due to the 
small values ctab = 0.8oaa and ctbb = 0.88ctaa- The lo- 
cal minima depend on the attraction between B particles 
within the A particles' FCS. Their maximum distance 
must be rs oab + caa = 1-8oaa> which is slightly less 
than 2oaa- Neverthelss, Figs. [5jb) and^b) show that 
this sensitivity to the attraction between pairs of B parti- 
cles within the FCS of the A particle has no consequence 
for their dynamics. 



IV. DISCUSSION 

We have demonstrated that structure and dynamics of 
LJ systems in the dense fluid phase are well reproduced 
by introducing a shifted-forces cutoff [l?} at the first min- 
imum of the radial distribution function (a distance that 
approximately equals the radius of the sphere for which 
the density is equal to the overall density in the uniform 
fluid). This means that for LJ liquids - and presumably 
also for other simple liquids - including all interactions 
from the particles within the FCS is all that matters for 
getting the correct physics. The forces from particles 
beyond the FCS only change the structure marginally, 
and their contribution to the thermodynamics can there- 
fore be taken into account with high accuracy by simple 
mean-field corrections. 

It is important to emphasize that this paper focused 
on the uniform condensed phase, far from the gas phase 
and the critical point. The dominant role of the forces 
within the FCS can only be expected to apply for the 
condensed phase, where a coordination shell is well es- 
tablished. When the density is decreased, one gradually 
enters a region of phase space of more gas-like behavior; 
here the cluster picture of the virial expansion gradually 
becomes more relevant and useful. For a simple LJ liquid 
our simulations suggest that at typical temperatures this 
transition takes place when pa 3 ~ 0.6. 

With the demonstration of the importance of the FCS 
for simple condensed fluids one may ask to which extent 
this property applies also for complex molecular liquids 
such as water? Some indication that this may be the case 
is given in Ref. H^, where the free energy of a (TIP3P) 
water molecule is partitioned into chemical associations 
with proximal inner-shell water molecules and classical 
electrostatic/dispersion interactions with the remaining 
outer-shell water molecules. The calculated free energy 
is in excellent agreement with free energy per particle 
of (TIP3P) water, but it should be noted that of course 
not all interactions beyond the FCS are ignored in this 
approach. 

The traditional discussion of the roles of attractive ver- 
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sus repulsive forces [1, H, Q, which was recently rein- 
vigorated by Berthier and Tarjus [H, [l6[ , does not ask 
the right question. According to the presently proposed 
"FCS focus" the reason that the traditional separation 
into repulsive and attractive interactions fails (severely 
for the viscous dynamics, but actually also for the stan- 
dard LJ liquid) is not per se that attractions are ignored. 
Rather, the dynamics and (first-order) perturbation ex- 
pansion fail because at ordinary densities ignoring the at- 
tractions fails to take into account all interactions from 
particles within the FCS - at extremely high densities 
the WCA approximation works well because here it does 
take into account all forces from particles within the FCS. 

More generally, the physical basis for the traditional 
perturbation theory 0, [t| is called into question by our 
findings. The idea behind perturbation theory is that re- 
pulsive and attractive forces give rise to separate, clearly 
identifiable contribution to the equation of state and to 
the free energy. Perturbation theory for the free energy 
is based on the physical picture that, roughly speaking, 
the repulsive forces decrease the entropy whereas the at- 
tractive forces decrease the energy. This picture captures 
the general qualitative behavior of many condensed fluids 
[37j . We find, however, that in order to obtain an accu- 
rate particle distribution all forces from particles within 
the FCS must be included, independent of whether these 
forces are repulsive or attractive. As a consequence, we 
suggest that perturbation theory will be much more use- 
ful if the reference state as a minimum takes into account 
all interactions from the FCS. 

The mean-field term is the first term in the Zwanzig 



high-temperature expansion [fjj. It is exact if there is 
a well-defined excluded volume in the high-temperature 
limit. This is, however, not the case. Our simulations at 
the high-temperature state point (3) demonstrate that 
even in this case a necessary and sufficient condition for 
obtaining the correct dynamics and thermodynamics is 
to include all forces from the FCS - and only these. The 
size of the FCS is density dependent, and no well-defined 
"high-temperature limit" exists. The Zwanzig reference 
system is a hypothetical reference system like the ideal 
gas state. As is well known, such systems are very useful 
for obtaining a qualitative understanding of the liquid 
state, as well as quantitative information, as has been 
demonstrated in numerous cases. Understanding the lim- 
itations of reference systems is nevertheless important 
for obtaining an accurate physical understanding of the 
physics of simple, condensed fluids. 
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